PURE POINT SPECTRUM OF THE FLOQUET HAMILTONIAN 
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Abstract. We prove that the 1 — d quantum harmonic oscihator is stable under spatially 
localized, time quasi-periodic perturbations on a set of Diophantine frequencies of pos- 
itive measure. This proves a conjecture raised by Enss-Veselic in their 1983 paper [EV] 
in the general quasi-periodic setting. The motivation of the present paper also comes 
from construction of quasi-periodic solutions for the corresponding nonlinear equation. 
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1. Introduction and statement of the theorem 



The stability of the quantum harmonic oscillator is a long standing problem since 
the establishment of quantum mechanics. The Schrodinger equation for the harmonic 
oscillator in (in appropriate coordinates) is the following: 

i=l * 

where we assume 

for the moment. We start from the 1 dimensional case, n = 1. (1.1) then reduces to 



The Schrodinger operator 



J?=(-;^+-^) (1-3) 



is the 1-d harmonic oscillator. Since H is independent of t, it is amenable to a spectral 
analysis. It is well known that H has pure point spectrum with eigenvalues 

An = 2n+1, n = 0, 1.., (1.4) 

and eigenfunctions (the Hermite functions) 

where Hn{x) is the n*^ Hermite polynomial, relative to the weight {Hq{x) — 1) 
and 

/oo 
e-""^ Hm{x)Hn{x)dx 
-oo (1.6) 

Using (1.4-1.6), the normalized solutions to (1.1) are all of the form 

oo 

ijix,t) = Y,(^nhnix)e'^' (^|a„|2 = l), (1.7) 

n=0 

corresponding to the initial condition 

oo 

V'(x,0) = ^a,/i„(x) (J]|a,|2 = l). (1.8) 

n=0 

The functions in (1.7) are almost-periodic (in fact periodic here) in time with frequen- 
cies A„/47r, n = 0, 1... 
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Equation (1.2) generates a unitary propagator U{t, s) = U{t — s,0) on L^(]R). Since 
the spectrum of H is pure point, Vtt e -£'^(K), Ve, 3i?, such that 

inf||C/(t,0)«||i,2(|,|<^)>(l-e)||«|| (1.9) 

by using eigenfunction (Hermite function) expansions. 

The harmonic osciUator (1.3) is an integrable system. The above results are classical. 
It is natural to ask how much of the above picture remains under perturbation, when 
the system is no longer integrable. In this paper, we investigate stability of the 1-d 
harmonic oscillator under time quasi-periodic, spatially localized perturbations. To 
simplify the exposition, we study the following "model" equation: 

d 1 " 
-^Ol'^ = + ^^)'^ + ^\^o{x)\^^cos{ujkt + 4>k)ilJ, (1.10) 

k=l 

on Ci(R, L2(R)), where 

< 5 < 1, 

^ = K}fc=i e [0,27r)^ 
<^={'^fc}Lle[0,27r)^ 

ho{x) = ' is the 0th Hermite function. 



(1.11) 



In particular, we shall study the validity of (1.9) for solutions to (1.10), when U is the 
propagator for (1.10). The method used here can be generalized to treat the equation 



where V is Cq^ in x and analytic, quasi-periodic in t. 

The perturbation term, 0{S) term in (1.10) is motivated by the nonlinear equation: 

-i^^iP=^{-^+x^)iP + MiP + 5\ipfiP (0<5«1), (1.12) 

where M is a Hermite multiplier, i.e., in the Hermite function basis, 

M = diag(Mn), MnGR, 

oo 

Mu = Y^ Mn{hn, u)hn, for all u e L^{R). 

n=0 

Specifically, (1-10) is motivated by the construction of time quasi-periodic solutions to 
(1.12) for appropriate initial conditions such as 

i> 

V'(x,0) = ^Cfc,/ife,(x). (1.13) 



In (1.10), for computational simplicity, we take the spatial dependence to be |/io(ic)P 
as it already captures the essence of the perturbation in view of (1.12, 1.13, 1.5). 
The various computations and the Theorem extend immediately to more general finite 
combinations of hk{x). 

The Floquet Hamiltonian and formulation of stability 

It follows from [Y2, 3] that (1.10) generates a unique unitary propagator U{t,s), 
t, s e M on L^(M), so that for every s e M and 

uoeB' = {feL\R)\\\f\\l,^ J2 < oo}, 

\a+f3\<2 

u{-) = U{-,s)uo e C\R, L2(M))nC°(M, B^) (1.14) 

is a unique solution of (1.10) in L^(R) satisfying u{s) = uq. 

When I' = l, (1.10) is time periodic with period T = 27r/a;. The 1-period propagator 
U{T + s, s) is called the Floquet operator. The long time behavior of the solutions 
to (1.10) can be characterized by means of the spectral properties of U{T + s, s) 
[EV,Ho,YK]. Furthermore the nature of the spectrum of U is the same (apart from 
multiplicity) as that of the Floquet Hamiltonian K [Yl] : 

d 1 

^ " '"^90 ^ 2^~dx^ ^ "^^^^ ^ '^l^o(a^)P COS0 

on L^(R) (8) -£'^(T), where L^(T) is I/^[0, 27r) with periodic boundary conditions. 

Decompose L^(]R) into the pure point Tipp and continuous He spectral subspaces of 
the Floquet operator U{T + s, s): 

L^{R) =Hpp®Hc. 

We have the following equivalence relations [EV,YK]: u e Hpp{U{T + s, s)) if and only 
if Ve > 0, 3i? > 0, such that 

^T^h\\U{t,s)u\\Lmx\<R) > (i-e)l|w||; 
and u e V.c{U{T + s, s)) if and only if Vi? > 0, 



1 

limt^ioo- t^^'l|t^(*',s)w||i2(|^|<^) = 0. 



(Needless to say, the above statements hold for general time periodic Schrodinger 
equations.) 

When z/ > 2, (1.10) is time quasi-periodic. The above constructions extend for 
small 5, cf. [Be, E, JL] leading to the Floquet Hamiltonian K: 

fc=i fc=i 
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on L^{R) L^iT"), cf. [BWl]. This is related to the so called reducibility of skew 
product flows in dynamical systems, cf. [E]. We note that the Hermite- Fourier func- 
tions: 

e-'''-'^hj{x), neZ", cj) e T'^, j G {0, 1...} (1.16) 

provide a basis for L'^{R) ^ L^iT"). 

We say that the harmonic oscillator H is stable if K has pure point spectrum. 
Let s G M. This implies (by expansion using eigenfunctions of K) that given any 
u e L2(M), Ve >0,3R> 0, such that 

inft||C/(t, s)m||l2(|^I<^) > {1- e)\\u\\,a.e.(f), (1.17) 

cf. [BWl, JL]. So (1.9) remains valid and we have dynamical stability. 
We now state the main results pertaining to (1.10). 

Theorem. There exists Sq > 0, such that for all < S < Sq, there exists Q C [0, 2tt)'^ 
of positive measure, asymptotically full measure: 

mes Q — > {27r)'^ as 5 — > 0, 

such that for all a; G O, the Floquet Hamiltonian K defined in (1.15) has pure point 
spectrum: cr{K) — Upp. Moreover the Fourier- Hermite coefficients of the eigenfunc- 
tions of K have subexponential decay. 

As an immediate consequence, we have 

Corollary. Assume that Q is as in the Theorem. Let s G M. For all u> & Q, all 
u G L^(]R), all e > 0, there exists R> 0, such that 

inft||t/(t,s,0)tt||L2(|a;|<ii) > (1 - e)||tt||, a.e.0, (1.18) 

where U is the unitary propagator for (1.10). 

We note that this good set f2 of a; is a subset of Diophantine frequencies. This is 
typical for KAM type of persistence theorem. Stability under time quasi-periodic per- 
turbations as in (1.10) is, generally speaking a precursor for stability under nonlinear 
perturbation as in (1.12) (cf. [BWl, 2]), where M plays the role of ui and varies the 
tangential frequencies. The above Theorem resolves the Enss-Veselic conjecture dated 
from their 1983 paper [EV] in a general quasi-periodic setting. 

A sketch of the proof of the Theorem. 

Instead of working with K defined on L^(]R) ®L'^(T^) directly, it is more convenient 
to work with its unitary equivalent H on P{1J^ x {0, l---}), using the Hermite-Fourier 
basis in (1.16). We have 

1 5 

i? = diag (n-a; + j + -) + -W^® A (1.19) 



on £^(Z'^ X {0, 1...}), where W acts on the j indices, j = 0, 1, 2..., 




VJ+7 
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_e 20+i') forj+/>l; 



(1.20) 



A acts on the n indices, n G Z*^, 



Ann' = 1, \n- n\ei = 1, 
= 0, otherwise. 



(1.21) 



The reduction from (1.15) to (1.19-1.21) is performed in section 2. The main work 
is to compute W, which involves integrals of products of Hermite functions. We 
will explain shortly this computation, which is independent from the main thread of 
construction. 

The principal new feature here is that W is long range. The j^^ row has width 
0{^/j) about the diagonal element Wjj. It is not and cannot be approximated by a 
convolution matrix. The potential x'^ breaks translational invariance. The annihilation 
and creation operators of the harmonic oscillator a= ■^(^+^); ct* — ^) ■> 

satisfying [a, a*] = 1, are generators of the Heisenberg group. So (1.19) presents a new 
class of problems distinct from that considered in [Bl-3, BWl, 2, EK, Kul, KP]. 

The proof of pure point spectrum of H is via proving pointwise decay as — ?/| — > oo 
of the finite volume Green's functions: {H\ — E)~^(x, y), where A are finite subsets of 
X {0, 1...} and A x {0, 1...}. We need decay of the Green's functions at all 

scales, as assuming E an eigenvalue, a priori we do not have information on the center 
and support of its eigenfunction ip. The regions A where {Ha — E)~^ has pointwise 
decay is precisely where we establish later that ip is small there. 

For the initial scales, the estimates on Ga{E) = {Ha — E)~^ are obtained by direct 
perturbation theory in 5 for < 5 <^ 1. For subsequent scales, the proof is a multiscale 
induction process using the resolvent equation. Assume we have estimates on Ga' for 
cubes A' at scale L' and A is a cube at a larger scale L, L ^ L' . Intuitively, if we 
could establish that for most of A' C A, Ga'{E) has pointwise decay, then assuming 
we have some a priori estimates on Ga{E), we should be able to prove that Ga{E) 
also has pointwise decay. 

There are "two" directions in the problem, the higher harmonics direction n and 
the spatial direction j. The off-diagonal part of H is Toeplitz in the n direction, 
corresponding to the discrete Laplacian A. Since the frequency uj is in general a 
vector (if > 2), n ■ uj does not necessarily — > oo as |n| — > cxo. So the n direction 
is non perturbative. We use estimates on Ga' and semi- algebraic techniques as in 
[BGS,BW1] to control the number of resonant A', where Ga' is large, in A. 

In the j direction, we do analysis, i.e., perturbation theory. This is the new feature. 
From (1.19) and Schur's lemma, A|| = 0{1). So the £^ norm of the perturbation 

does not decay (relative to eigenvalue spacing) in j. However when 5 = 0, if is 
diagonal with eigenvalues n ■ uj -\- j and eigenfunctions 5n,ji the canonical basis for 



e{Z^ X {0, 1...}). We have 



\\[W®A]5nA = 0{^) (j>l), 
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which decays in j. 



This is intuitively reasonable, as W stems from a spatially localized perturbation 
from (1.10). As j increases, The Hermite functions hj become more extended, c.f. 
(1.5). So the effect of the spatial perturbation should decrease as j increases. 

Assuming uj is Diophantine: 

\\n ■ <^||t > -j — (c > 0, n 7^ 0, q: > 21/), 

where || • ||t is the distance to the nearest integer, this enables us to preserve local 

eigenvalue spacing for A which are appropriately proportioned in n, j. This in turn 
leads to decay of Green's functions. Combining the estimates in the n and j directions, 
we obtain estimates on the Green's function at the larger scale L. 

Integrals of products of Hermite functions. 

From (1.15, 1.19), computation of W involves computing the following integrals: 



/. 



oo 
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hQ{x)hm,ix)hn{x)dx 
=== e ^''^H^{x)Hm{x)Hr,{x)dx, m, n = 0, 1..., 



where Hm, are respectively the m*^, n*^ Hermite polynomial, Hq{x) — 1. 
Let 



/oo 
e-''^' Hl{x)Hm{x)Hn{x)dx. 
-oo 



-oo 

The idea is to view Hq{x) as Ho{\/2x), i.e., the 0**^ Hermite function relative 
to the weight e~^^ and to use the generating function of Hermite polynomials to 
reexpress 



Hm{x)Hn{x) = aiHi{V2x). (1.23) 



£=0 

We then have 



I = ao J [Ho{V2x)fe-^''^dx 

using (1.6). This computation is carried out in section 2, (2.7-2.10), recovering an 
apparently classical result, which could be found in e.g., [GR, PBM]. 

More generally, we are interested in computing 

/oo 
e-'^^'^ Hp{x)Hq{x)Hm{x)Hn{x)dx, p, q, m, n = 0, 1... (1-24) 
-oo 



which are needed for the nonhnear equation or if we consider more general perturba- 
tions of the harmonic osciUator. FoUowing the same hne of arguments, we decompose 
Hp{x)Hq{x) into 

P+Q 

Hp{x)Hg{x) = J2 biHt{V2x). (1.25) 
e=o 

Combining (1.23) with (1.25), and assuming (without loss of generality), p+q- < m+n, 
we then have 

p+q 
e=o 

where q = JZ,[He{V2x)]^e-^^" dx. 

The computation for general p, q is technically more involved and is carried out in 
[W]. Unlike the special case p = g = 0, we did not find the corresponding result for 
general p, q in existing literature. 

The computation of / in (1.24) is exact (see (2.10)), reflecting the integrable nature 
of the quantum harmonic oscillator. The proof of the theorem is, however, general. 
It is applicable as soon as the kernel W satisfies (1.20). Following the precedent 
discussion on / for general p, q, and using properties of the Hermite series (cf. [T] and 
references therein), one should be able to extend the Theorem to V which are Cq^Iu x 
and analytic quasi-periodic in t, leading to perturbation kernels in the Hermite-Fouricr 
basis satisfying conditions similar to (1.20) in the j direction and exponential decay 
condition in the n direction. 

When the perturbation V is independent of time and is a 0*^ order symbol, satisfying 

\d''V\<Cail + \x\)-°', a = 0,1... (1.26) 

the corresponding Schrodinger equation has been studied in e.g., [BBL, KRY, Z], where 
it was shown that certain properties of the harmonic oscillator equation extend to the 
perturbed equation. The spectral property needed for the construction here is more 
detailed and stringent. Hence it is reasonable to believe that the set of potentials V 
will be more restrictive than that in (1.26) 

Some perspectives on the Theorem 

The Theorem shows that for small 6, there is a subset O C [0, 2-71")'^ of Diophantine 
frequencies of positive measure, such that ii u> E then (1.9) is satisfied. Hence 
spatially localized solutions remain localized for all time. It is natural to ask what 
happens if the forcing frequencies u are in the complement set, co & Q'^. 

If uj is rational, the perturbing potential V is bounded and has sufficiently fast decay 
at infinity, it is known from general compactness argument [EV] that the Floquet 
Hamiltonian has pure point spectrum. In our example, this can be seen as follows. In 
(1.10) restricting to periodic perturbation {u = 1), it is easy to see that for 

Va;, A = {noj + j + z)~^W A, where $52; = 1 is compact. (1-27) 

Assume to is rational: uj = p/q, {q 0). Since Hq = nuj+j has pure point spectrum 
(with infinite degeneracy) and the spacing between different eigenvalues is 1/q, (1.27) 
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implies that H has pure point spectrum. When oj is irrational, Hq typically has dense 
spectrum. No conclusion can be drawn from (1.27). 

It is worth remarking that (1-27) holds for all scalar u). In the quasi-periodic case, 
w is a vector, the compactness argument breaks down. The proof of Lemma 3.5 in the 
present paper is a replacement. 

If V is unbounded, wc have a different situation. The results in [HLS, GY] combined 
with [YK] show that for the following unbonded time periodically perturbed harmonic 
oscillator: 

duL 1 

i— = + x^)u + 2e{siivt)xiu + ijy{t,x)u, x = {xi, ...Xn) e^"" (1.28) 
where F(t, x) is a real valued smooth function of (t, x), satisfying 

V{t^2T^,x) = V{t), \V{t,x)\<\x\ asa;^oo, \d'^V{t,x)\<Coc, |ci;| > 1, 

the solutions diffuse to infinity as t — > oo. More precisely, for all uq G L^(M'^)ni7^(M"), 
for any i? > 0, the solution ut satisfies 

1 f^ 

hm -/ cit||Wt||L2(|a;|<ii) = 0. (1.29) 
T->-±oo I Jo 



In (1.28), V = \ (periodic), w = 1 , a; G O'^, (1.29) is an opposite of (1.9). However 
the perturbation is unbounded. Moreover the proof in [GY] uses in an essential way 
that the potential is linear at infinity, hence positivity of the commutator: [^^, x{\ = 1. 

In the exactly solvable case where the time periodic perturbations is quadratic in 
the spatial coordinates, it is known that the Floquet Hamiltonian exhibits a transition 
between pure point and continuous spectrum as the frequency is varied [Col]. The 
perturbation there is again unbounded. 

Some related results. 

To our knowledge, when a; G O is nonresonant, there were no results in the literature 
on the perturbed harmonic oscillator equation of type (1.10), even in the time periodic 
case, i.e., u G [0,27r). The main difficulties encountered by the traditional KAM 
method seem to be (i) the eigenvalue spacing for the unperturbed operator docs not 
grow, Afe+i — Afc = 1, (ii) the perturbation W in the Hermite basis has slow decay 
(1.20). 

When the eigenvalue spacing for the unperturbed operator grows: \Ej^i — Ej\ > 
{(3 > 0), which corresponds to a potential growing faster than quadratically at 
infinity, and when the perturbation is periodic in time, related stability results were 
proven in [DS]. In [Co2], under time periodic perturbation and replacing W in (1.20) 
by a faster decaying kernel, hence decaying norm in j, which no longer corresponds 
to the physical case of harmonic oscillator under time periodic, spatially localized 
perturbation, stability results were also proven. Both papers used some modified 
KAM method. 
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Motivation for studying (1-10) 

As mentioned earlier, the motivation for analyzing (1.10) partly comes from the non- 
linear equation (1.12). In [Bl-3, EK], time quasi-periodic solutions were constructed 
for the nonlinear Schrodinger equation in with Dirichlet or periodic boundary con- 
dition 

^-^^= (-A + M)V' + (^|V'P^V', (peN+;0<5<l) (1.30) 
ot 

where M is a Fourier multiplier; see [Kul, KP] for the Dirichlet case in M with a 
potential in place of M. In [BW2], time quasi-periodic solutions were constructed for 
the nonlinear random Schrodinger equation in 

d 

i—^|;=(-eA + V)^/; + 6\^{;\^Pu, (p e N+ ; e, < (5 < 1), (1.31) 

where V = {vjljgzd is a family of random variables. 

The proofs in [Bl-3, BW2] use operator method, which traces its origin to the 
study of Anderson localization [FS]. This method was first applied in the context of 
Hamiltonian PDE in [CW]. The proofs in [EK, Kul, KP] use KAM type of method. 

In (1.30) (specializing to 1-d), the eigenvalues of the linear operator are n^, so 
En+i — En ^ n, the eigenfunctions e*"'^, however, are extended: |e*"^| = 1 for all x. 
Let us call this case A, where there is eigenvalue separation. In (1.31), the eigenvalues 
of the linear operator form a dense set, the eigenfunctions, on the other hand are not 
only localized but localized about different points in Z*^ from Anderson localization 
theory, see e.g., [GB, GK]. This is case S, where there is eigenfunction separation. 
The existence of time quasi-periodic solutions, i.e., KAM type of solutions in ^ is a 
consequence of eigenvalue separation; while in B, eigenfunction separation. 

Equation (1.10) and its nonlinear counterpart 

-i^^^ = + ^')V' + MV^ + (^IV'I'^V, (p e N+; < 5 « 1), (1.32) 

where M is a Hermite multiplier, stand apart from both (1.30, 1.31). It is neither 

A, nor B. There is eigenvalue spacing, but it is a constant: A^+i — = 1. In 
particular, it does not grow with n. The eigenfunctions (Hermite functions) h^, are 
"localized" about the origin. But they become more extended as n increases because 
of the presence of the Hermite polynomials, cf. (1.5). This in turn leads to the long 
range kernel W in (1.20) and long range nonlinearity in (1.32) in the Hermite function 
basis, cf. [W]. 

From the KAM perspective a la Kuksin, this is a borderline case, where Theorem 
1.1 in [Ku2] does not apply. The more recent KAM type of theorem in [EK] does not 
apply either, because W is long range and not close to a Toeplitz matrix (cf.(1.20)) for 
the reasons stated earlier. These are the features which make (1.10, 1.32) interesting 
from a mathematics point of view, aside from its apparent relevance to physics. 
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2. The Floquet Hamiltonian in the Hermite-Fourier basis 



Recall from section 1, the Floquet Hamiltonian 

d 1 ^ 
^ = ^ E '^'^ + 2 ^~ + + (^l^ola;)!' (2.1) 

on L2(M) ® L2(T^), where 

< 5< 1, 

Wfc e [0,27r), k = 1, 

c/>fc e [0,27r), /c = l,..,z. (2-2) 

As mentioned in section 1, ho{x) is the 0*^ Hermite function, O*'^ eigenfunction of the 
1-d harmonic oscillator and more generally, 

An = 2n + 1, n = 0, 1... (2.3) 
hn{x) = ^^e-^'/^ n = 0, 1..., 

where Hn{x) is the n**^ Hermite polynomial, relative to the weight e~^^ (Hq(x) = 1) 
and 

00 

2 



e ^ Hm{x)Hn{x)dx 

(2.4) 

=2"n!x/7r5r„n 

Integral of products of Hermite functions 

We express (2.1) in the Hermite function basis for small 5 and compute the integral: 



/oo 
hl{x)hm{x)hn{x)dx 
-00 

1 f°° 2 

/ e'^"^ H^(x)Hm(x)HJx)dx, m, n = 0, 1..., 



(2.5) 



V2^ 

The integral of the more general product 



00 



e Hp{x)Hq{x)Hm{x)Hn{x)dx 

(2.6) 

= (i^p(a;)i^q(a;)e"^^ Hm{x)Hn{x)e~''^), p, q, m, n = 0, 1... 
is done in [W]. 
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We use generating functions of Hermite polynomials to find of (1.23) as follows. 
Since 



oo 



n=0 

oo 



^2sx-a'' 



V —Hm{x), 

ml 



(2.7) 
(2.8) 



m=0 



which can be found in any mathematics handbook (cf. [CFKS, T] for connections with 
the Mehler formula), multiplying (2.7, 2.8), we obtain 



^-^ n\m\ 

n,m 

oo / t+s \£ oo 

p=0 



(2.9) 



From (1.23), we are only interested in the coefficient in front of Ho{^/2x). So we set 
£ = 0. To obtain ao, we equate the coefficient in front of f^s"^. Comparing the LHS 
with RHS of (2.9), n, m must have the same parity, otherwise it is 0. We deduce 



(_1)^ 

an = — T ■ (n + m)l, n, m same parity 

2^(m±ny^ ^ ' ^ ^ (2.10) 

= otherwise, 

by taking 2p = n + m, which is the only contributing term. 

Taking into account the normalization factors in the third equation of (2.3), we then 
obtain 



Lemma 2.1. 



Wmn = / hl{x)hm{x)hn{x)dx 



— oo 



(m + n)! [tt 
2 



Let 



2-+-V^ (^)! 
=0 



m, n same parity 
otherwise. 



(2.11) 



n — m 



assuming n > m, without loss. When N ^ 1, 



IWrr 



< 



1 ^-fcV2Ar 



■Nl 



■s/2N{N + k)\{N -ky. 



N 



(2.12) 

(2.13) 
(2.14) 
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Proof. We only need to obtain the asymptotics in (2.13, 2.14). This is an exercise in 
Stirhng's formula: 

..= (^)"v^(l + ^ + ^ + ...) (2.15) 

or its log version 

logn! = (nH — ) logn — n + log \/27r + ... (2-16) 
Here it is more convenient to use the latter. Using (2.12, 2.16), 



(m + n)\ , (27V)! , 
loe — ■ — - — = loe — loe: 2 

=iVlogiV - iV + - log2 + 0(iV-^). 



So 



(m + n)! 



using (2.17). Hence 

(-l)TiV! 



m 



l + O 



^y2N{N + k)\{N -k)\ 



A^> 1 



(2.17) 



(2.18) 



which is (2.13). Using the fact that 

n\ = \/27rn ^- j e'*'" 

with 

i2;rT:i<^"<iL'^°^"^^"^'' 

and applying the inequalities (with x = k/N): 

(l){x)'^= (l+x) log(l + x) + {1- x) log(l - x) >x'^ 

for all x e [0, 1) and (f){x) > ax^ with a > 1 for x G [7/10, 1), we obtain (2.14). (When 
X = k/N = 1, (2.14) follows by a direct computation using Stirling's formula.) □ 

From (2.14), the matrix element W^n has subexponential decay 

\Wmn\<-^e-^-^ (iV=^»l) (2.19) 



when \m — n\ > y/N. When \m — n\ < y/N, we only have the estimate 

\Wmn\ < (2.20) 

V N 



Hence is a matrix with a slowly enlarging region of size 0{VN) around the principal 
diagonal Wnn, the P norm of W is of C(l), but the local £^ to norm is of order 
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0(1/ VN). These new features will need to be taken into account when we do the 
analysis in sections 3, 4. 

(2.1) in the Hermite-Fourier basis 

In the Hermite-Fourier basis, e-*""^/ij (a;), n e Z"", e T^, j e {0, 1...}, the Floquet 
Hamiltonian, which is unitarily equivalent to the K defined in (1.15) is then 

1 S 

H = (iiag{n-Lo + j+ -) + -W^A (2.21) 

on ^{^" X {0, 1...}), where W is the matrix operator defined in (2.11), acting on the 
J indices, j = 0, 1, 2..., A acts on the n indices, n G Z'^, 

A„„/ = 1, In — ril\n = 1, 

(2.22) 

= 0, otherwise. 

Let H = H — 1/2 and rename H\ let 5 = 5/2 and rename 5, 5. We then have 

r\pf 

H = diag(n ■ u; + j) + SW ^ A (2.23) 
on £2(Z^ X {0, 1...}), with W, A defined in (2.11, 2.22). 

3. Exponential decay of Green's functions at fixed E: estimates in 9 

Let H be the operator defined in (2.23), i.e., 

H = diag{n-u; + j) + 6W^A (3.1) 
on e^iZ" X {0, 1,...}), where 

n e Z^, j e {0, 1,...}, < 5< 1. (3.2) 
W acts on the j indices, j = 0, 1, 2..., 

w (-1)"^ 0"+/)! 

Vr,,/ = . , ., — - ■ — r— T7 — , 7, 7 same parity, 

2-?+-? VTlf! (^)! (^33) 

=0 otherwise. 



Write 



When J > 1, 



J = ^T— , k = ^r— , (assume j > j ). 



\Wjr\ < 



+ k)\{J -k)\ 



1 fc^ .X 

< ^e-- (3.4) 
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from (2.14). A acts on the n indices, n e Z'^, 



Ann' = 1, \n-n'\ii = 1, 

(o.o) 

= 0, otherwise. 

In view of the Theorem, our aim is to prove that on a good set of a;, H has pure 
point spectrum. To achieve that goal, we add a parameter 9 {9 eM.) to H: 



H{9) = H + 9 

= diag {n-io + 9 + j) + 6W®A. 



(3.6) 



We consider a sequence of finite volume Green's functions at fixed E: 

Ga{9,E) = {H^{9)-E)-\ (3.7) 

where A are cubes in Z'^ x {0, 1...}, k / v. {0, 1...} in an appropriate way, Ha{9) 
are H{9) restricted to A. 

In this section, i? G M, a; G [0, 27r)'^ arc fixed and we assume that c<; is a Diophantine 
frequency. We do estimates in 9. Specifically, we prove (inductively) that for any A 
large enough, away from a set of 9 of small measure in M, ||Ga(^, -E")!! is bounded 
and \Gp^{9,E){x, y)\ has subexponential decay for \x — y\ ~ linear scale of A. (For 
precise statement, see Proposition 3.10.) As mentioned in section 1, the proof is a 
combination of a non perturbative part in the n direction and a perturbatuve part in 
the j direction. We note that this set of bad 9 depends on E, u), the estimate on the 
measure is, howevere, uniform in E and u for Diophantine uj. 

In the next section (section 4), we eliminate the E dependence by excluding double 
resonances and converting the estimate in 9 into estimates in uj in the process. The 
conversion is possible because lv, 9 appear in (3.6) in the form n ■ co + 9. 

Below we start the induction process. To simplify notations, we extend i7 to a 
linear operator on ^^(Z^+i) with Wjj> as in (3.3) for j, f G {0, 1...} and Wjj> = 
otherwise. 

3.1 The initial estimate (0*^ step). 
For any subset A C Z'^'^^, we define 

Hj,{9){n,j;n',j') = H{9){n,j;n',j'), {n,j;n',f) G A x A, 
= 0, otherwise. 

Let Aq = [— J, J]'^'^^ for some J > to be determined. For a fixed E, we study the 
Green's function Gao{9,E) = {H\q{9) — E)~^ by doing perturbation theory in 5. We 
have 



Lemma 3.1. Assume < (5 <^ 1. For any fixed a , Q < a < 1/4, there exists J G N 
such that the following statement is satisfied. Let Aq = [—J, J^'^^ . There exists a set 
B{Ao,E) in R, with 

mes B{Ao,E) <e--^^^\ (3.9) 
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such that if 9 e M\B(Ao,-E), then 

||GA,(^,E)||<e^^ (3.10) 
\G^,{e,E){n,r,n'3')\ < e-K-'^)-(-''^')l^'\ |(n, j) - (n',/)| > J/10. 

(3.11) 

Proof. Since this is the initial estimate, we do perturbation theory in 5. Let k > 0. 
Let B(Ao, E) be the set such that if ^ e ^(Ao, E), then 

\n-uj + j + e - E\<2k (3.12) 

for some (n, j) G Aq. Clearly 

mes B(Ao, -E) < 4|Ao|k 

= 4(2J+l)^+i«. ^^-^^^ 

Let 

i?A„(^)='diag(n-c^+j + ^), (n,j)eAo (3.14) 

be the unperturbed diagonal operator. Since \\5W ® A||^2 < 0(5), if ^ ^ B(Ao,-E'), 
then 

\\GK,{e,E)\\ = \\{D^,{e)-E + SW^A)-^ 

/ _i (3-15) 

for K 5. 

From the resolvent equation 

Gao{0, E){n,j; n'f) ^ {n ■ u + j + 9 - E)-'[{SW ® A)Gao{0, E)]{n,r, n'f) (3.16) 
for {n,j) 7^ {n',j'). Hence 

\GA,{e,E){nJ;n'j')\ < 0(l)5«-^ if \{nj) - (n',j')| > J/10 (3.17) 

Let 

J = I log^l (hence 6 = e'^), (3.18) 
K = e--'\ (3.19) 

(3.13-3.17) then imply (3.9-3.11). □ 

3.2 A Wegner estimate in 9 for all scales. 

We now state an apriori estimate in 6 for \\{H/^{6) —E)~^\\ valid for all finite subsets 
A C I/^'^ and all 5. This estimate will be useful in the induction process. Following 
the Anderson localization tradition, we call it a Wegner estimate. 
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Proposition 3.2. For any E & W, and any finite subset A in 1^'^^ , the following 
estimate is satisfied for all k> 0: 

mes {e\ dist {E, a{HA{e))) < k} < 2\A\k. (3.20) 



Proof. Let A/., k = 1, |A| be eigenvalues of Ha{9 = 0). Then 

{^1 dist {E,a{HA{e)))<K} = u^,^j^^{e\\E-e-Xk\<K}. (3.21) 

It follows that the measure of the left side is bounded by 2|A|«;. □ 

Our goal now is to obtain inductively the equivalent of estimates (3.9-3.11) for larger 
subsets A, Z'^^^ D A D Aq. We note that in proving Lemma 3.1, we did perturbation 
theory about the diagonal operator 

Dao (0) = diag {n-u; + j + d)\ ^^j^eAo (3-22) 

using the smallness of S. This was sufficient for one initial scale. For subsequent scales, 
however, we need more detailed information on the spectrum of Ha- 

3.3 Local spectral property of Ha- 

We assume that u is Diophantine, i.e., 3 c > 0, a > 2i/, such that 

\n-uj + j\>-^ (3.23) 

for all n G Z^\{0}, all j G Z. In this subsection, we make statements which hold for 
any fixed 6. We look at finite subsets A C Z^+^, such that j ^ 0, if {n,j) £ A. When 
j ^ 1, this is the perturbative region. 

Proposition 3.3. Assume oj satisfies (3.23). Let < /3' < l/Sa and 3/4 < (3 < I. 
Let A be a rectangle centered at (N, 2L) e Z"^ x Z, where iy is identified with iy x {0}; 

A = (TV, 2L) + \-L^\l^'Y X [-L^, L^\ C Z'^+\ (3.24) 

Assume 

\L\ > 1. (3.25) 

For any fixed 6* G M, the eigenvalues Xn,j{0) of Ha{0) satisfy 

\KAO) - K',j'm > (n, j) ^ (n'Z); (3.26) 

the eigenfunctions (f)n,j may be chosen such that 

\\(l>n,j-Sn,j\\< (3.27) 
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Remark. It is crucial to note that the estimates (3.26, 3.27) are independent of $ and 
the specific co satisfying (3.23). In Lemma 3.5, we exploit further the consequences of 
(3.26). 

Proof of Proposition 3.3. Let (n, j), (n', j') G A, (n, j) ^ [n' , j'). Then for all 9, the 
difference of the diagonal elements 

\n-u+j+e-in'-u+j'+e)\ = \in-n')-uj+ij-j')\ > > for < /?' < l/5a, L » 1, 

(3.28) 

from (3.23). Use as approximate eigenfunctions Snj with approximate eigenvalues 
Xn,j{0) = n ■ uj + j + 9, {n, j) e A, and let 



{Ha - ~Xn,j)Sn,j = V'- (3.29) 



Then from (3.7), 



ip{n\j') = SWj'j, \n-n'\ = 1, (n, /) G A 
= 0, otherwise. 



(3.30) 



and 



from (3.4). 



= C»(|L|-i/^) (3.31) 



Equations (3.28-3.31) imply that Xn,j{9) = n-uj+j + d is an approximate eigenvalue 
of H\ to 0{\L\~^/'^). This can be seen as follows. Let A = Xnj{9). Assume 

dist {X,a{HA{e)))>0{\L\-y^). 

Take any /, ||/||^2(a) = 1, / = Y^Cmi^rn, E |c^P = 1, where -i/^m are eigenfunctions of 
Ha. Then 

{Ha - A)/ = ^{Em - X)Cm'lpm, 

where is the corresponding eigenvalue for i/j^^. Hence 



\\{Ha - A)/|| = ^JJ2(^rn - A)2|c^|2 > 0{\L\-'/^). (3.32) 

This is a contradiction if / — ^n,j • 

Hence (3.28) gives (3.26). (3.27) follows from (3.26, 3.31) and standard perturbation 
theory, see e.g., [Ka]. □ 



3.4. The first iteration (1®* step). 

We now increase the scale from J to J*-^, where C > 1 (independent of S) is the 
geometric expansion factor, which will be specified in section 3.5. Recall that J is the 
initial scale, large enough so that (3.18) holds. Hence on Aq = [—J, JY~^^, (3.9-3.11) 
hold. 
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Let 

All? [-J^^J^Y+K (3.33) 

Our aim is to prove the analogue of (3.9-3.11) when Aq is replaced by A. The general 
strategy in going from scale J to scale is to distinguish the region near j = 0, 
where we use the estimate from the previous scale, here (3.9-3.11) and non perturba- 
tive arguments and the region away from j = 0, where we use Proposition 3.3 and 
perturbation theory. The general iteration strategy here is similar to that in [BWl] . 

Toward that end, we define 

T""^ {(n,j)eA|U|<2J-l}. (3.34) 

Let 

Ao(n, j) = (n,j) + [-J, J]^+\ (n,j) e T (3.35) 
be cubes of the previous scale. Let 

A,(n,2j) = (n,2j)+[-J^',/yx[-/,/] (0 < < l/5a, 3/4 < /3 < 1), (n, 2j) G A\T, 

(3.36) 

be cubes of type (3.24). We cover A with Aq, A* cubes, i.e., T with Aq, A\T with A*. 
Ga are then obtained by using resolvent equation, (3.9-3.11), Propositions 3.2 and 3.3. 
Wc implement this strategy in our first iteration. This iteration is special, as in the 
tube region T, wc use smallness of 5, cf. (3.18). 

We need the following notion of pairwise disjointness. Let Sk, k = 1, be finite 
sets , Sk 7^ Sk' ii k ^ k' . Let S — {Sk}^^i- If Sk H 5"^/ ^ 0, V/c 7^ A;', then we say 
there is 1 pairwise disjoint set in S. More generally, if 3/i, I2, Ip, Ip n Ipi — 0, if 
p 7^ p\ {Ip}p^i = {1, 2, K} such that Sk n Sk' 7^ if and only if A;, k' e Ip for some 
p. Then we say there are P pairwise disjoint sets in S. 

Lemma 3.4. Let Fq he a covering of T with Aq cubes defined in (3.35). Assume uj 
is Diophantine satisfying (3.23). Fix E , a {Q < a < 1/ A) as in Lemma 3.1. For all 6, 
there exists at most 1 pairwise disjoint A e Fq, such that 

dist{E,a{Hj^{e)))<e-^\ (3.37) 
Moreover if Aq E Fq and Aq n A = 0, 

dist{E,a{HAM))>^- (3.38) 

Proof. Let Aq G Fq. If \Xn - E\ < e~-^^ for a G a{HAo{0), then since \\5W (g) A|| < 
C6 = Ce-J, 

\n ■ uj + j + e - E\ < e'-^" + Ce--^ < 2e~-^'' 
for some (n, j) G Aq. Thus if both Aq and A satisfy (3.37), then 

\{n - n') ■ + j - j'\ <\n ■ + j + - E\ + \n' ■ + j' + e - E\ <4e~-^" (3.39) 
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for some (n, j) G Aq and {n' e A. This implies Aq n A 7^ for large J, since 
otherwise the left side is larger than c/ {2vJ'^)°^ by the Diophantine condition (3.23), 
which is a contradiction if C <C J*^/ log J. 

If Aq e To is disjoint from A e Fq, then for any (n, j) e Aq and for some {n',j') e A, 
\n-uj+j + e-E\ > \{n-n')-uj+j-j'\-\n'-uj+j' + e-E\ > ^^Jjcy -2e~-^". (3.40) 

Since \\5W ® A|| < C5 = Ce"-^, (3.38) follows. □ 

If Aq n A = 0, then (3.10, 3.11) are available. Let be a covering of A\T. For 
A* e r*, we need 

Lemma 3.5. Let < < 1/5q; and 3/4 < /3 < 1. Let A* he the rectangle: 

KM) = (n,2j) + [-f\fr X [-/,/], M) e A\r. (3.41) 

For a fixed E E M., there exists W G M, with mes W < C|A*|e~"^'' such that for 
9 e R\W, J > 1, 

|[^A.W-^]-'Kj>",/')l<e"^^'"'"""'+^^, 

i\n' - n"\ > /'/lO, or \j' - j"\ > f/10) (3.42) 

ll[^A.(^)-^]-'ll<e-^''''- (3.43) 



Proof. For any fixed 9, E, there exists at most 1 bad site b = {n',j') € A*, such that 

1 

|n ■cc;+/ + 6'-E| < 



This is because if there were {n' {n" e A*, {n' ^ {n'\j") such that 

1 



\n' -uj + j' + 9-E\< 



\n"-u+j" + 9-E\< ^ 



then 

|(n'-n")- a; +(/-/)!< ^ 



^•1/5' 



which contradicts the Diophantine condition (3.23) on ut: 



- -Ik ^ ^ » (3.44) 



for j > J ^ 1, using (3.41). 
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To obtain (3.42, 3.43), we first consider H\^\i,{9) and make estimates on [Hj^^\i,{6) — 
E]~^{n' ,j';n" We then use Proposition 3.2 and the resolvent equation to obtain 
(3.42). We prove below that 

\\[HKAb{0) - E\-%.^,. < 2//5^ (3.45) 

for J 3> 1. To obtain estimates on the matrix elements, we use weighted £^ space and 
show that [H^^\p{9) — E]~'^ remain bounded. (3.45) is then the special case when the 
weight equals to 1. 

Toward that end, for any p G Z^'^^, we define \p\y — |(pi, 0)| and \p\i = 

1(0, ...,0,p^^i)\. Let p — {pa}aei''+^ be a family of weights, such that 

p^(p) = e'^""'''^^ (aeZ'^+i), VpeZ^+\ (3.46) 
where j is as in (3.41). Hence Wp, q G Z^'^'^, 



Pa {P)Pa{q) = e • e 



< e 



P ll-^^ •1/2 



(3.47) 



for all Pa & p- 

Let be the diagonal part in (3.8). To arrive at (3.42), we consider the deformed 
operator HA,\biO): 

HAAbm?, q) = [p-'HAM^)pa]{p, q) 

=Dp, + Sp-\p)pa{q){W ® A)p„ (3.48) 

=Dpq + S(W A)pq 

where p, q & A*\6 and 

{W A)p/= p-\p)Pa{q){W ® A)p,. (3.49) 

To prove boundedness of [H^Abi^) ~ we use resolvent series and perturb 

about the diagonal. We have formally 

[HaM^) - ^]"' =(D -E)-^+ 5{D - E)-\W 0A){D- E)-^ 

+6^{D - E)-\W ® A){D - E)-\W A){D - E)'^ (3.50) 
+... 

Let 

yV'^= {W A){D - E)-\ 



Then 



W{n',j'; n",j") - Wj>j>>A^,^.{n" ■ u + j" + 6 - E)-\ \n' - n"| = 1 
= otherwise. 
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(3.51) 



We bound W using Schur's lemma [Ka] : 

(sup J2 mn',j';n",j")\y/\sni^ Yl IMn' ,j';n" (3.52) 



< 



where (n',/), {n",j")eK\h. 

J2 mn',j';n",f)\= J] J] \WjijiiAr,in"{n" ■ u; + f + 9 - E)-\ 

n",j" \n"-n'\ = l |i"-2i|<j/5 

(n',/) e K\h. From (3.47, 3.4) 

,200 (3-53) 
where we used j' + /' < Aj + 2j'^ < 9j /2 for /, j" e A*\6. So 



J] |W(n',/;n",/')l<^ E E |n" • + /' + ^ - 

n",j" "^-^ \n"-n'\ = l\3"-2j\<3f> 



(3.54) 



<^(//^ + logj)- 

Let aj/''== n" ■ u + j" + 6 ^ E. To arrive at (3.54), we used the fact that |aj//| = 
\n" ■ uj + j" + 9 - E\ > for all {n",j") e A*\6 and that ap - aq = p - q, for all 

p, q. 

Using (3.53), we have 

J2mn',j';n",j")\<0{l)j'/\ (3.55) 



Substituting (3.54, 3.55) into (3.52), we have 



i^/s 1/9 



So the Neumann series in (3.50) is norm convergent: 

WIHaM^) - <\\{D- E)-mi + 6\m + 6'\\Wr + ...) 
<2\\iD-E)-'\\ 
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which is equivalent to 

\\Pa{HAM^) - Ey'Pa'W < 2//^ (3.56) 
For each pair p, q E A*\6, we can then always choose a so that 



[HAM^)-E]-\n',j';n",j")\ < C//5e-[l-'-"l+^]. (3.57) 



Write Go for [Hj^^\jj{9) — E] ^, G for [Ha,{0) — E] ^. From the resolvent equation: 

G = Go + GoH^Gq + GoHhGHhGo, 

where Hj, = H^^ — HA,\b- The matrix element Hi,{p^ q) = 0, unless p = b or q = b. 

Using this, (3.57), the Wegner estimate (3.20) on G with k = e"-^" '"' (0 < < l/5a) 
and self-adjointness, we obtain (3.42, 3.43) □ 



We now write the estimate at scale J^. Assume 5, J satisfying (3.18), so that 
Lemma 3.1 holds. Let Ji = J'^ , C > 1, the same geometric expansion factor as before. 
Let A = [-Ji, Jil^+K We have 

Lemma 3.6. Assume uj is Diophantine satisfying (3.23) and < 5 <^ 1 is the same 
as in Lemma 3.1. For any fixed cr, < u < 1/ba, there exists B(A, E) in with 

mes B{A, E) < e'-^^^^ (3.58) 

such that ifde R\B{A,E), then 

\\GA{e,E)\\<e''? (3.59) 
\GA{e,E){n,j;n'f)\ < e-K-'^)-(-''^')l^'', 

for all {n,j), (n'j') such that \ {n,j) — {n',j')\ > Ji/10, 

(3.60) 

provided the expansion factor C satisfies 1 < C < (3' /a, < /3' < l/5a is as in (3.41)- 

So we have the same estimate as in Lemma 3.1 at the larger scale Ji = J'-^ {C > 1). 

Proof. This is similar to the proof of Lemma 2.4 in [BWl]. So we summarize the main 
steps. We prove (3.59, 3.60) using the resolvent equation and cover A with cubes of 
types Aq, a* defined in (3.35, 3.36), i.e., T, defined in (3.34), with a covering Fq of 
Aq's and A\T, a covering F* of A^'s. From Lemma 3.4, for all fixed E, all 9, there 
exists at most 1 pairwise disjoint Aq G Fq on which (3.10, 3.11) do not hold. We use 
Proposition 3.2 on this Aq, (3.10, 3.11) on all other Aq. 

For a given A*, let Wa, be the set such that (3.42, 3.43) hold if 6* G M\VFa,. Let 
W = UWa^ , where the union is over all possible A* with centers in A\T. 

mesVv<0{l)Ji e <e 1 ■ Ji , 
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where the first term is an upper bound on the number of possible A* with centers in 
A\r multiphed by the volume of A*. For 9 e we use (3.42, 3.43) on A*. 

Using the resolvent equation, combining (3.10, 3.11, 3.42, 3.43) and (3.20) with 
K = e~'^^ (0 < cr < I/Sck) on the only bad Aq, we obtain (3.59, 3.60). (For more 
general iterations using the resolvent equation, cf. proof of Lemma 3.8, in particular 
(3.72).) Combining the measure estimate from (3.20) with k — e~"^i and the above 
measure estimate on W, we obtain (3.58), provided C < (3' /a. □ 

For the induction process to follow, it is convenient to define the following. For any 
fixed CT, < cr < 1/5q; and any given box A C 'L^'^^ of side length 2J + 1, we say 
Ga{0, E) at fixed (6*, E) is good if 

||GA(^,i?)||<e^'' 

|GA(^,£^)(n,j;ny)|<e-l(-'^)-('^''^')l''', 

V(n, j), (n'/)such that |(n, j) - (n',/)| > J/10. (3.61) 

Otherwise, it is had. We also define 

>V(A)^^*U^A., (3.62) 

where the union is over all possible A^^ of the form (3.36) with centers in A\T and T 
is as defined in (3.34), so that for 9 E R\W(A), (3.42, 3.43) are vahd for all A* with 
centers in A\T. We say A* is good, if (3.42, 3.43) hold. Otherwise A* is bad. 

3.5 A large deviation estimate in 9 for the Green's functions at fixed E at 
all scales. 

We now increase the scale from Ji to Ji (1 < C < (3' /cr^ the geometric expansion 
factor will be determined here). Our task is again to derive estimates (3.58-3.60) for the 
cube [- Jf , Jf ]^+^ starting from the estimates (3.58-3.60) for the cube [-Ji, Ji]^+^ 

For simplicity of notation, we rename Ji, J and [—Ji,JiY~^'^, A in this section. 
As in the first iteration, we distinguish the tube region T, defined as in (3.34) with the 
new J. We cover T with Aq's defined in (3.35) with the new J, and A\T, A* defined 
in (3.36). We note that Aq arc at scale J with centers in T, while A* are at scales 
from Jf^' to J*^^ (0 < < 1/5q;, 3/4 < /3 < 1) with centers away from T. 

As in the first iteration, we use the resolvent equation to obtain estimates on Ga 
from estimates on Ga„ and Ga, . Let Fq be a covering of T and F* of A\T. For 
9 e M\>V(A), (3.42, 3.43) are valid on all A* G F*. So for any fixed 9 e ]R\W(A), we 
only need to control the number of pairwise disjoint bad Aq boxes on which estimate 
(3.61) is not available. In particular, we need the number of such bad boxes to be 
<^ J^, the linear scale of the box A. (This is intuitively clear, as otherwise without 
further detail on the location of the bad boxes, we could not accumulate decay at the 
linear scale as in (3.60).) Recall that for the first iteration, there is at most 1 such 
(pairwise disjoint) bad box. 
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Lemma 3.7. Assume u is Diophantine, satisfying (3.23) and Lemma 3.6 is valid on 
cubes 

Ao(0, j) = (0, j) + [-J, JY+\ Vj e [-(2 J - 1), 2J - 1]. (3.63) 
Then for all fixed 9 

#{(n, j) G T|Ao(n, j) zs a 6a(i &oa; } 
<0(l)j5(.+i) ^ (jC)i- <^ jc^ (3-64) 

6?/ choosing 5(iv + 1) < C < J""/^ (t) < cr < l/Saj. 

Proo/. Write Aq for Ao(0,j). We first replace the estimate \\G^q{9,E)\\p^p < e-^^ in 
(3.59) by the estimate on the Hilbert-Schmidt norm: 

||G'Ao(^,i?)||HS<e^r (3.59') 

This leaves the measure estimate in (3.58) unchanged for Ji ^ 1 (cf. proof of Lemma 
3.6). 

Define 

A"^' U B{Ao{0,j)). (3.65) 

je[-{2j-i),2j-i] 

Since the conditions on the Green's function in (3.59', 3.60) can be rewritten as poly- 
nomial inequalities in 6 by using Cramer's rule, A is semi-algebraic of total degree less 
than 

(2 J + 1)2(^^+1) . (2 J + 1)2(^^+1) . (4J + 1) 

< 0^(1) J^(^+^), (3.66) 

where the first factor is an upper bound of the degree of polynomial for each entry 
of the matrix Gaq{9,E), the second is an upperbound on the # of entries of each 
G\q{9,E) plus the one for the Hilbert-Schmidt norm, the third is the # of diff'erent 
matrices Gaq's. For more details, cf. the proof of Lemma 2.6 in [BWl]. A is therefore 
the union of at most 0^{1)J^^'^~^^^ intervals in M by using Theorem 1 in [Ba] (see also 
[BGS], where the special case we need is restated as Theorem 7.3). 

For any fixed 9 eR, let 

I = {ne[-J^,J^]''\n-uj + 9 eA}. (3.67) 

Then Ao(n, j) is a bad box if and only if n e / and, therefore, n is in one of the 
intervals of A. But each interval does not contain two such n if a; satisfies (3.23) and 
\I\ < Oi/(l) J^^'^"'""^-' by virtue of (3.58). This is because for Diophantine ui satisfying 
(3.23), if there exist n, n' e [-J^, J^]", n n', then 

\in -n')-u;\> » e-^^''^ (3.68) 

Hence each interval can contain at most 1 integer point in [— J*^, J*^]^. We therefore 
obtain (3.64). □ 

For any fixed 9 E M\W(A), W(A) defined as in (3.62), the only bad boxes are of 
type Aq. Lemma 3.7 shows that there are only few (of order ((J*^)^~) bad Aq boxes 
in A. The following iteration lemma will enable us to obtain estimates (3.58-3.60) for 
G\ at scale J*-^. 
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Lemma 3.8. Fix b & (0, 1/8) and assume r satisfies 3/4 + b < t < 1 — b. Suppose 
M, N are integers satisfying 

<M < 2N^. (3.69) 

Let A = [— A^, A^]'^"'"^. Assume for a// A C A with diameter L, the Green's function 
Gx{E) = (i?^ — E)~^ at energy E satisfies 

\\Gk\\ < (3.69') 

Let A' be cubes of side length 2M. We say that A' is good if in addition to (3.69'), 
the Green's function exhibits off-diagonal decay: 

\GA'{E){x,y)\<e-\^-y\''" 

for all x,y e A' satisfying \x — y\ > M/10. Otherwise A' is bad. Assume for any 
family T of pair-wise disjoint bad A' cubes in A, 

#JP < N\ 

Under these assumptions, one has 

\Gj,{E){x,y)\<e-\--y\''" 
for all x,y & A satisfying \x — y\ > N/10, provided N is sufficiently large, i.e., N > 

No{b,T). 

Proof. The proof is similar to the proof of Lemma 2.4 in [BGS]. As we will see from 
(3.73, 3.74), because of the conditions on 6, r, it only needs a one step iteration. 

To estimate Gx{E){x,y), x, y G A, \x — y\ > A^/10, let Q be cubes of side length 
4M, we make an exhaustion {Si{x)Yi=o of A of width 2M centered at x as follows: 

S-^{x)''^ 0, 
So{x) A'(x) n A, 
Siix)-"^ U^es,_,g(y)nA, 

for 1 < i < £, where £ is maximal such that Si{x) ^ A. 

We say an annulus Ai — Si{x)\Si-i{x) is good if AiHJ-' — 0. Let Ai{x), Ai+i{x), Ai. 
be adjacent good annuli and define 

U = Ul+J,A,{x). 

Let d^S-i{x) = {x} and 

d,Sj{x) = {ye Sj{x)\3z e A\Sj{x), \y - z\ = 1} 

for j > 0. 

26 



By construction 



dist {d^Si-i,d^Si+s) = 2M{s + 1) (s > 0). 
For any subset B C A, let Hb be defined as in (3.8). Let 

Tb *== Ha — {Hb © Ha\b)- 
Assume x & B, from the resolvent equation, 

GA{E){x,y) = GB{E){x,y) 

+ Gb{E){x,z)Tb{z,z')Ga{E){z' ,y). (3.70) 

zee 
«'eA\B 

The proof follows the same line of arguments as in the proof of Lemma 2.4 in [BGS] 
by iterating (3.70). There are two modifications: 

• In view of F^, with matrix elements Fb(z, z') — 5VF(j, j')A(n, n'), z = {n,j) e B, 
z' = {n',j') e A\B, where A{n,n') as defined in (2.22) and 

\W{j,j')\ < e 20+j') < e 4iv— 
from (2.14) for K. j + j' < 2N, we define 

dSl'\x) = {ye Si{x)\dist{y,d.Si{x)) < N^^/^""}, 
dSPix) = {ye A\Si{x)\dist{y,d.Si{x)) < N^^/^""}. 

We note that for y e dsl^\x) U dsl'^\x), dist{y,d^S^{x)) < Ar"/i6 ^ m, the size 
of A' cubes, and for j, j' such that \j - j'\ > N^^/^^, 

|W^(i,/)l«e-^'''. (3.71) 

• For all X E with dist(x, d^Si-i) > M/4, there exists x' E U such that A'{x') C 
U and dist(a:,a*A'(x')) > M/5. We estimate Gu{E){x,y) with x e dsf\, y e 
dS\^2^ using A\x') C U. 

To obtain subexponential decay of off-diagonal elements GA{E){x,y) we proceed 
as in [BGS]. This entails to estimate iteratively Gs^^.{E){x, z), where z G dSn-j^ and 
Am = Sn,\Sni-i is a good annulus. Let Ui < nii < Ui+i, so that all the annulus in 
Smi\Sni are bad and all the annulus in Sni_^.l\Smi — U are good. 
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Using the resolvent equation to relate Gg^.^^ z)-,^ ^ dSn}+i with Gs^. {E){x, y), 

y e dSn} , we have 

Gs^^^^iE)ix,z) = Yl GuiE)iz,w)Tuiw,w')Gs^^^^iE)iw',x) 



well 



J2 Gu{E)iz, w)Tuiw, w')Gs^^^^ {E){w' , x) + o(e-^''') 



J2 Gu{E){z, w)ru{w, w')Gs^^ {E){y, x)rs^^ {y, y')Gs^^^^ {E){y', w') 

(3.72) 

where we used (3.71, 2.22). This is the analogue of (2.25, 2.26) in [BGS]. Iterating 
(3.72) and taking the log log, lead to the conditions 

loglog[(e^'')^/^] < logloge^^'\ (3.73) 

logloget^'^'W^-^')] > logloge^''', (3.74) 

for N 1, where (3.73) originates from the estimates on ||Gs„.^J| and (3.74) from 
the decay estimates on Gu- (3.73, 3.74) in turn lead to 3/4 + 6 < r < 1 - 6 for M, iV 
satisfying (3.69) (cf.[BGS]). □ 

In order to apply the above lemma, we need to convert the covering of A = 
[-J^, J'^]'^+i which is of diverging scales from Jf^' (0 < /5' < l/5a) to J^^ (3/4 < 
(3 < 1) to a covering of a single scale J^^. Let A' be cubes at scale J^^ (3/4 < < 1), 

^{n,j) = {n,3) + [-J<'f',J^^Y+\ {n,j)eZ^+\ (3.75) 
Let W be as in (3.62). We define 

A' bad (for a fixed 9 e R\W(A)), if A' n A 7^ 0, where A is a bad Aq box, i.e., (3.61) is 
violated. (3.76) 

Let r' be a covering of A with A' boxes defined in (3.75). Lemma 3.7 gives 

Lemma 3.9. Fix a G (0, l/5a) as in Lemma 3.6. For any fixed 9 G M\W(A), V has 
at most 0(1) J^^^^+i) = {.J^Y~ < (^{^ + 1) < C < J<^/2; pairwise disjoint bad A' 
boxes. On the good A' box, we have 

||G^,(^)||<e(^^^)^ 
|GA.(^;n,j;ny)|<e-l("'^>("''^'')l''', 

V(n, j), {n'j') such that \{n,j) - {n',j')\ > J^^/10, (3.77) 
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provided 

max(5(i/ + 1), < C < ^ < r'''. 

zpa a 

So (3.61) is satisfied for the good A'. 

Proof. From (3.64), there are at most 0{1),J"'^^^^^ points (n, j) G T, such that Ao(n, j) 
is a bad box, where the 0{1) only depends on the dimension p. Since the A' boxes 
are at scale J^^ » J, the scale of Aq boxes, there are at most 2^+^ non intersecting 
A' boxes intersecting a given bad Aq box. Hence for any fixed 9 G M\W(A), V has at 
most (9(1) J^(^+-'^^ pairwise disjoint bad A' boxes with a slightly larger 0{1). 

On the good A' boxes, we use the resolvent equation to obtain (3.77) as follows. 
We first estimate ||Ga'(^)||- Let e > be arbitrary. We have 

Ga{E + ie){x,y) = Gu(x)iE + ie){x,y) 

+ Yl Gu(x){E + ie){x,z)ru(a:){z,z')GA{E + ie){z',y), 
zeu{x) 

z'€A\U{x) 

where Tu(^^) is as defined above (3.70), U{x) is either A* or Aq and x G U{x). Write 
U for U{x). For y' G dM = {y" e U\3z G A\U, \z -y"\ = 1}, x and y' satisfy the 
distance condition in (3.42) iiU — and the condition in (3.11) iiU — Aq. It is easy 
to see that for all a; G A, such a U (x) exists. 

Summing over y, we have 

J2\GAiE + ie)ix,y)\< ^ \\Gu\\ 
y y€U{x) 



+ 2^+V^(^+i) 6-^5 sup ^ \Ga{E + te){w, y)l 



where we used (2.14, 2.22, 3.10, 3.11, 3.42, 3.43) and the conditions a, (3' < l/5a, 
a > 2v. Taking the supremum over x and [/(x), we obtain the bound on ||Ga'(6')|| in 
(3.77), provided C > P'/2Pa, J > 1. 

Applying the resolvent equation one more time, similar to the proof of Lemma 3.8, 
and using the bound on ||Ga'(6')||, we obtain the off-diagonal decay in (3.77). □ 

Let < /?' < l/5a, (a as in (3.23)), 3/4 < /? < 1. Choose a G (0, l/Sa) satisfying 
+ 1) < P' /o-, which leads to 

0<a< ^ < 



5(zy + l) 25a(iy + l)' 

Choosing appropriate C large enough so that both Lemme 3.8 and 3.9 are available, 
we then arrive at the following estimate for Ga{6) = {Ha{0) — E)~^ valid for all 
A = [—J, JY'^^, with J large enough and any fixed E. 
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Proposition 3.10. Assume uj is Diophantine, satisfying (3.23). For any < P' < 

l/5a, fix 

0<a< -. 

20(z/+ 1) 

Then for all < 5 1, there exists Jq such that the following statement is satisfied 
for all A = [- J, J]''+^, with J > Jq: There exists B{A, E) in R, with 

mesB{K,E)<e-^'''\ (3.78) 
such that ifOe R\B{A,E), then 

\\GA{e,E)\\<e'\ (3.79) 

|GA(^,£;)(n,j;ny)|<e-l(-'^-)-("''^')l''', 

for all {n,j), (n'j') such that \{n,j) — {n',j')\ > J/10. 

(3.80) 



Proof Choose C > 40(zv + 1) such that 

|<C.</*'<1. 

Assume < 5 <^ 1. Let Jq = | log . For the scales Jq to Jq , we use perturbation 
theory in 6 as in Lemma 3.1, with a further lowering of 6 if necessary in order that 
(3.9-3.11) hold for all scales in [Jq, Jq]. For the scales > Jq , Lemme 3.8 and 3.9 are 
available, in view of the choice of C. 

Assume (3.78-3.80) hold at some scale J » 1. Using Lemma 3.9 in Lemma 3.8, 
we obtain the corresponding estimate at scale J^. Hence the proposition holds by 
induction. □ 



4. Frequency estimates and the elimination of E. 

We now convert the estimate in 9 in (3.78-3.80) for fixed E, uj into estimates in uj 
for fixed 9 (6* = 0), eliminating the E dependence in the process by excluding double 
resonances. This is the key estimate leading to the proof of the Theorem in section 5. 
We will first define a double resonant set in (a;, ^) e x R. 

We use 2 scales N, iV_with N < N, loglogiV < logiV. Let AnU) = [-N, N]''+^ + 
(0,j), j e Z, AjY = [—N,N]'^'^^. Assimic A^ is resonant at ^ = for some E, i.e., 
E is close to some eigenvalue of H/^^ (see (4.1) below). We prove in Lemma 4.1 that 
with further reductions in the frequency set, the AAr(j) boxes with \ j\ < N and are at 
appropriate distances from A^ are all nonresonant (at ^ = 0). 

Toward that end, let DC{N) be the Diophantine condition to order N, i.e., if 
LU e DC{N), then 

||n • a;||T > ^ (« > Vn e [-N, iV]^\{0}, 
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where as before || • ||t is the distance to the nearest integer. Denote H/^^ by Hff. We 
define V{N, iV) c x R as 

p(iV, N)= y {(a;, 6*) e T'^ X M\3E, such that 

{jez:|Aiv(j)nAiv#0} 
||(i/^(a;,0)-£;)"'|| >e^' and ^ G -B(Aiv(j), (4.1) 

and for fixed < a < 1 and C > 1, 

S{N)= U {V{N,N)(\DC{N)xR)}, (4.2) 

where N x A^*-" means aA^*-^ < -/V < N'^ /a. We note that a is kept fixed in this section. 

The set S{N) is the double resonant set in (cu, 9) restricted to w G DC{N), S{N) C 
X M. It is the projection of the set S{N), S{N) C T'^ x M x M, i.e., S{N) = 
Proj T-xR('S'(iV)) with 

S{N) = [j {V{N, N) f]{DC{N) xRx R)}, (4.3) 

where 

V{N,N) = {iu;,d,E) e X R X R| ||(i?^(a;, 0) - £;)"^|| > e^' 

and 9 e B(Aiv(j), E) for some AnU) n Ajv ^ 0}, (4.4) 

Lemma 4.1. Let N be sufficiently large, C > 1 and S > sufficiently small. Let 
Q,N be the union of D,^ over N x N*-^ , where C DC{N) is the set of to, such that 

• there exists E & such that 

||(iy^(a;,0)-E)-'||>e^'; (4.5) 

• there exist j e Z, Aiv(j) n Ajv 7^ 0, and I e Z", \£\ x exp[(fog A/")^] = K such that 

|(iyA,(,-)-E)-^(a;,^a;)(m,m')| >e-l™— 'I''' (4.6) 

for some m, m' E An (j) , \m — m' \ > N/10. 
Then the set Qn satisfies 

mes < exp[-^(fogiV)2]. (4.7) 



The proof of Lemma 4.1 uses the foUowing decomposition lemma [B3, Lemma 9.9]. 
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Lemma 4.2. Let S C [0, l]^'^ be a semi- algebraic set of degree B and mes2i/<S < 
rjjlogB <^ logl/?7. Denote by [x^y] e [0,1]^ x [0,1]^ the product variable. Fix e > 
^i/2iy_ r/ien there is a decomposition 

with Si satisfying 

\Proj^Si\<B^e (4.8) 
and S2 satisfying the transversality property 

mes^(52 n L) < B^e-^T]^/^", (4.9) 
for any v- dimensional hyperplane L such that 

ma^^ \ProjL{ej)\ < (4.10) 

where ej are the basis vectors for the x- coordinates. 

For our usage, the variable {x, y) G in the lemma will be (cj, 9) after identifying 
9 with (^,0) e W. 

Proof of Lemma 4-1- Let {u!,9) G S{N), then 9 G B(^AN{j), Ek) for some Aiv(j) and 
some eigenvalue Ek of Hj^ satisfying \Ek—E\ < e~^ . This is because the perturbation 
e~^ essentially preserves the condition of the definition of a bad set, see (3.8-3.10) 
(c.f., proof of Lemma 4.1 [BWl]). By the Diophantine restrictions on the frequencies 

mes S{N) < C>(l)e-^^''iV^("+') (4.11) 

from Proposition 3.10, where the third factor in the RHS is an upper bound on the 
number of AAr(j) and the number of eigenvalues of Hfj. 

Since S{N) is semi-algebraic with total degree at most N'-"^^ (Ci > 1), cf., proof 
of Lemma 3.7, S{N) its projection onto T*^ x M is also semi-algebraic with degree at 
most N'^^" (C2 > Ci) [Chap 9, B3]. (4.8) and Lemma 4.2 then conclude the proof by 
taking e x exp[-(logiV)2] (cf., also [BW2]). 

Here we also used the fact that the sum of the measure estimate in (4.9) over all 
L corrsponding to the hyperplanes (c<;,£ ■ u) with \i\ x exp[(log A?")^] is much smaller 
than exp[— (log A'")^] for the corresponding B, rj and e > 0. □ 

Lemma 4.1 shows that the boxes Aiv(-^,j) = + [-N,N]'^+^, where \£\ x 

exp[(logiV)^], j G [-A^'^,A^'^] are non resonant with the box Apj- = [-iV, iV]''+^, 
N X N'-^ . Let T be the tube region, 

T={{n,j) G Z^+'\j G [-N^,N%. (4.12) 

The boxes Ajv(-^, j) have centers in T. 

We now exclude resonances of A^^ with boxes with centers in 7i^\T. Recall from 
section 3, that boxes with centers in Z^\T are in fact rectangles: 

A^n, 2j) = (n, 2j)+[-/',/']-x [-/,/] (0 < < l/5a, 3/4 < /? < 1), (n, 2j) G Z^T. 

(4.13) 

To exclude resonances of boxes A*(n, 2j) where max(|n |,|2j|) X exp[(logiV)2] with 
A^, we use direct perturbation. We have 
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Lemma 4.3. 

mes {u e [0,27r)''| dist {a{Hj^),a{HAj) < k, lj satisfies (3.23)} 

1 9 (4-14) 

<Cav|A;vI|A*|, (AC<exp[--(logiV)2]) 

o 

for any A* = A*(n, 2_;) in (4-13) satisfying \n\ > exp[i(log A^)^] and max(|n|, \2j\) x 
exp[(logA'')^], provided < f3' < min(l/20, l/Sa). 

Proof. Let A^^^, 0rn,fe be eigenvalues and eigenfunctions of H/^^: 
Write (f)m,k = S(m',fc')eA* (^m',k'Sm',k', with Yl |am',fe'P = 1- Then 

Am,fc = ^ |a^',fe'P(m' • a; + A;') + 5 ^ Clm',k'0.m",k"^m'm"Wk'k"- 

{m',k') (m',k'),(m",k") 

From (3.27) of Proposition 3.3 

J2 = 0(^). (4.15) 



The first order eigenvalue variation: 

d_ 

duo 



m'k'l'^'m ■ m' 



{m',k') 

= ^ \am',k'\'^rn- m+ ^ \am',k'\'^rn ■ {m' - m) 

{m',k') {m' ,k'),m'^m 

where we used (4.15) and the fact that \m — m'\ < 0{l)j^ in A*. Lowering P' to 
satisfy f3' <1 /20 if necessary, we obtain 

d 1 

m ■ g^^m,k > 2 > (4-16) 

for |m| > exp[i(logA^)2] > N'^ . 

Let iJ,£^i{u!) be the eigenvalues of i^jv- Then /x^,i(w) are piecewise holomorphic in 
each (jjp, p= l, and Lipshitz in u: \\iJ,e,i\\up < CN x A?"^. Using this and (4.16), 
we have that 

mes {a; e [0,27r)''| minmin \Xm,k{<^) - l^e,i{'^)\ < «} 

£,i m.k 



<C«|A^||A,|. 



□ 



Using Lemma 4.3, we arrive at 
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Lemma 4.4. Let N be sufficiently large, C > 2/(3' , < P' < min(l/20, 1/5q;) as 
in Lemma 4-3, S sufficiently small. Let be a subset of the Diophantine set defined 
in (3.23) with the properties: 

• there are N x N*-^ and E such that 

||(iyjv(c^,0)-£;)~'|| >e^'; (4.17) 

• there is a 

A(n,2j) = (n,2j) + [-/,/]-+!, (3/4 < /3 < 1) 
with \j\ > 12 and max(|n|, \2j\) x exp[(log A^)^] such that 

\{H^^r^,2j)-E)-\u;,0){m,m')\ > e"!— '1^'' (4.18) 

for some m,m' E A satisfying \m — m'\ > j^/10. 
Then 

mes n'j^ < e-^. (4.19) 

Proof. We cover A with A* of the form (4.13). We use a resolvent expansion similar to 
the proof of Lemma 3.9 to obtain the opposite of (4.18). For that purpose we need to 
take away an additional set of u> so that (3.42, 3.43) hold on all A* of the form (4.13) 
at ^ = 0. 

For any A*(n, 2j) such that |n| > exp[i(log A^)^], max(|n|, |2j|) x exp[(log A^)^], 
(4.14) is available. Take one such A* and let n = e~^^ » for N large. Assume 

(4.17) holds, so \E - Ek\ < < e~^^ for some eigenvalue E^ of i^jy-. Lemma 4.3 

then says that 

\\{E - Hj,J-'\\ < e'^^ll + 2e-^') < 2e^^ 
by taking away a set in uj of measure 

<C«|A^||A,| (4.20) 
Comparing (4.12) with (3.34) J = here. So 

||(£; - i^Aj-^ll < 2e^^ < e-^'^'"' = e^"''''^ 
provided C/3'/2 > 1 or C > 2//3'. Therefore (3.42, 3.43) hold on this A^. 

Multiplying (4.20) by the number of all possible A* (n, 2j) such that \n\ > exp[i (log A^)^] , 
max(|n|, \2j\) x exp[(log A^)^], we have that (3.42, 3.43) are available on all such A* 
after taking away an additional set in uj of measure <e~^. 

For A*(n, 2j) such that \n\ < exp[|(log AT)^], since max(|n|, |2j|) x exp[(log AT)^], 
\j\ X exp[(logAr)2]. So 

dist {E,a{HAj) X exp[(logA^)2], Vo; e [0,27r)^, 

ifE satisfies (4.17). This is because < 0{N) = <^ exp[(log A^)^]. 

The proof now proceeds as the proof of Lemma 3.9 and we arrive at the conclusion. 

□ 
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5. Proof of the Theorem 



We use Lemme 4.1, 4.4 to prove that the Floquet Hamiltonian K in (1.15), or 
rather its unitary equivalent H in (1.19) has pure point spectrum. Let VLn-, ^'n be the 
frequency sets as in Lemme 4.1, 4.4, N eN, sufficiently large. Define Q to be 

= T''\ IJ {flN U n'^), No > 1 (depending on 6), (5.1) 

Ar>iVo 

then mes Q — > {2n)'^ as 5 ^ 0. 

Fix u E fl. From the Schnol-Simon theorem [CFKS, S], to prove H in (1.19) has 
pure point spectrum, it suffices to prove that the generalized eigenfunctions have fast 
decay, hence are in More precisely, let be a non zero function on Z^^^^ satisfying 

{H-E)il; = 0, |V'(m)| < 1 + Iml'^o for all m e (5.2) 

where E is arbitrary and cq > is some constant. We will prove using Lemme 4.1, 4.4 
that ip has subexponential decay and hence is in £'^{Z'^'^^). 

We first verify that (4.5) is satisfied. This implies that (4.17) is also satisfied as 
they are the same condition. So we need to show that there is some box A^^ centered 
at 0, A;v = iV]^+i, for some TV x AT^ (C> 2/13', < (3' < min(l/20, 1/5q;) as in 
Lemme 4.1, 4.4), such that 

\\{Hf;,-E)~'\\ (5.3) 

For this we let 

r = {(n,j)eAAr|je[-iV^',iV^']}, 

C" > 1 to be determined from (5.8-5.11), C chosen to be > C' . Let Ajyc' be boxes of 
side length 2N'^' . We cover T with A^c> , An\T with A* of the form (4.13). 

From Proposition 3.10, Lemma 3.7, there are at most iV'^~'^° (tq > 0) pairwise 
disjoint bad A^c' boxes in T by taking C large enough. In Aj^\T, there is at most 1 
pairwise disjoint bad A* box by a further reduction in the Q set as follows. 

Assume 3(n, j) e A*, 3(m, k) e A'^, A* n A'^ = such that 
where Xnj, are eigenvalues of -f^A*, -f^A;- Prom Proposition 3.3 

So 

if n = m, since {n,j) (m, k). 
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We only need to look at the case n ^ m. This is similar to the proof of Lemma 4.3, 
except we look at the first order eigenvalue variation in the (n — m) direction. Let 

be the eigenfunction with eigenvalue Xn,j' 

HA^fpnJ — Xn,j4'n,j'i 

and 

(m',fe')eAi 

be the eigenfunction with eigenvalue ^: 

HA'^tpm^k = ^m,k'ipm,k- 

We have for the first order variation: 

d 

d 

= {n-m)-—[ ^ \an',j'\^n' - u - ^ \hm' ,k'\^m' ■ u] 

^ (n',i')eA. {m',k')eK 

d 

= (n — m) ■ — [(n — m) • a; + \an' {n' — n) ■ u 

J2 \bm',k'\^im' -m)-u;] (5-5) 

(m',fc')eA'^, m/m' 

= \n-mf + 0{j^) ■ f\n - m| + 0{^)k^'\n - m\ 

1| 12 

>-\n — m\ 
1 

where we used (3.27) of Proposition 3.3 and the definitions of A*, A'^ in (4.13). So 
mes {ijj\\\n,j{oj) — ^'m,k{'^)\ ^ i^i <^ Satisfies (3.23)} < Cn. 

Take k = e~(^°s Ar)^ ^^^^ denote the set Q'^ such that 3A* C Ajv, A'^ C A^, A*nA'^ = 
0, 3{n,j) e A*, 3{m,k) e A'^, such that (5.4) holds. Then 

mesn';,<e--^^'°^''^\ (5.6) 



provided 2//3' < C < logiV. 
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Let 

Q^'MqX [j Q%, (5.7) 

7V>iVo 

mes n (27r)'^ as 5 ^ 0. We now assume u E O,. Then in A^, N x ^ there are at 
most iV^""^" (0 < To < 1) bad Aj^c' boxes with centers in T and 1 bad A* box with 
center in Ajy\T. Since 

for C large enough, there has to be an annulus ^ at a distance x AT^ to the origin of 
thickness lOA^"-^ in T and of thickness /A in Ajy\T (as there is at most 1 bad A* 
box) devoid of bad points. 

Since {H - E)i/; = from (5.2), let 

(^A(c) - E)^ = (5-9) 

where A(c) is a box of type A = Aj^c' or A* centered in c. In view of (1.20, 5.2), the 
restriction on /? for A* in (3.36) 

if dist(m, ^A^-*^)) > \m — c|i, where 

aA(i) = {ye A\3z e Z'^\A, \z - y\i = 1} 

and 111 is as defined above (3.46). 

Let c e A, A(c) C A. So A(c) is a good box of type A^c' with centers in T or A* of 
the form (4.13) with centers in A\T . Estimates (3.42, 3.80) are available. From (5.9) 

^=(^A(c)-^)"'eA- 

Using (5.10), (3.42) or (3.80), we have that 

\m\<e~''\ (5.11) 

provided C > 2/(3' for aU c e A, such that 3A(c) C A. We now choose 2/(3' < C < 
C <S log A?", so that (5.8) is satisfied. (5.11) implies immediately that 

\^{m)\<e-^^ (5.12) 

for all m E A' (Z A, where A' is a smaller annulus contained in A of thickness 6A'"'-^ in 
T and thickness /5 in Ajv\T. 

Since dist A' to the origin x one can always find a square Aj^ centered at the 
origin, N x , such that 9Ajv C A', 

dist {dAj^, A'nT)> 2N^' , (5.13) 

dist {OAn, A' n T^) > A^^/10. (5.14) 
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Using (5.9) with Ajy- in place of A(c), we have 

Using (5.12-5.14), we have ||^jvll < e~^'^ by slightly increasing C if necessary. Since 
we may always assume ■0(0) = 1, (5.10, 5.11) give that 

\\{H^-E)-'\\>e''\ 



Subexponential decay of eigenf unctions of H. 

Let K = exp[(log N)"^]. Lemme 4.1, 4.4 then imply (via a proof similar to the proofs 
of Lemme 3.8 and 3.9) that the Green's function of the set 

C/='A2k(0)\Ak(0), 
where Al(0) = [— L,L]'^+-^, exhibits off-diagonal decay, i.e., 

|G'[7(E, m,n) I < e-l^-^l'^' m, n e t/, |m - n| > K/10. 

Let be defined as in (5.9) with U replacing A. Fix m E U. Assume dist (m, dU) > 
K/4:. For n such that dist {n,dU) > K/10, we use the bound in (5.10) with 
replacing Otherwise we use the fact that is polynomially bounded since ip is 
polynomially bounded. The equality 

iP{m) = ^G[/(£;,m,n)^(n), 

n 

then implies that 

|V'(m)| < e-l'"!" (0<K<l/4), 
provide K and thus |m| is large enough. □ 
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